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We investigate the time evolution of a system of quarks interacting with 
a and pion fields starting from an initial configuration consisting of a tube of 
hot quark plasma undergoing a boost-invariant longitudinal expansion. We 
work within the framework of the linear sigma model using classical transport 
equations for the quarks coupled to the mean-field equations for the meson 
fields. In certain cases we find strong amplifications of any initial pion fields. 
For large-radius tubes, starting from quark densities that are very close to 
critical, we find that a disoriented chiral condensate can form in the centre 
of the tube. Eventually the collapse of the tube drives this state back to 
the true vacuum. This process converts the disoriented condensate, domi- 
nated by long- wavelength pion modes, into a coherent excitation of the pion 
field that includes significant components with transverse momenta of around 
400 MeV. In contrast, for narrow tubes or larger initial temperatures, am- 
plification occurs only via the pion-laser-like mechanism found previously for 
spherical systems. In addition, we find that explicit chiral symmetry breaking 
significantly suppresses the formation of disoriented condensates. 
PACS numbers: 12.38.Mh, 11.30.Rd, 24.85.-hp, 25.75.-q 

I. INTRODUCTION 

Relativistic heavy- ion collisions provide a way to form regions of hot, dense hadronic 
matter. If the temperature is high enough, this matter is expected to be in a phase where 
chiral symmetry is restored and quarks are unconfined. Recently there has been much 
interest in the possibility that, as such a system cools, it could lead to regions in which the 
quark condensate is misaligned with respect to the physical vacuum. These regions can also 
be thought of as coherent excitations of the pion fields along particular directions in isospin 
space. Such a state is known as a disoriented chiral condensate (DCC) and a signal for its 
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formation would be anomalously large event-by-event fluctuations in the ratio of charged to 
neutral pions.[] 

Many of the studies of DCC's found in the literature use idealised geometries in order to 
simplify the calculations. These include uniform matter in a flnite box |Q, inflnite matter 
undergoing a boost-invariant expansion in one direction [0-0 and isotropically expanding 



inflnite matter ||T3|,|TJ,|12[ . An important exception is the work of Asakawa et al. who 
considered the classical evolution of the chiral flelds for cylindrical systems undergoing boost- 
invariant longitudinal expansion. 

In a previous paper (hereafter referred to as I), we examined the effects of flnite size 
on the evolution of the chiral flelds in the case of spherical systems. The framework used 
was the linear sigma model fl^ with explicit quark degrees of freedom. This model has 
also been studied by Csernai and Mishustin 0], for the case of inflnite matter expanding 
in one direction.^ Here we extend our studies to the case of cylindrical systems undergoing 
a boost-invariant longitudinal expansion. This geometry should be more relevant to ultra- 
relativistic heavy-ion collisions, and in particular to the central regions of such collisions 
where particle production is expected to be independent of rapidity [ pT[1 . 



As in Ref . , we assume a rapid quench that leaves the quarks out of thermal equilibrium 
with the chiral flelds. The subsequent evolution of the system is then described by the 
classical Euler-Lagrange equations for the flelds coupled to a relativistic transport equation 
for the quarks . We work in the classical limit where the transport equation reduces to a 
relativistic Vlasov equation for the distribution of the quarks in phase space. Further details 
of the model and method of solution can be found in I. 

For spherical droplets of hot quark matter, our studies in I showed that the quarks 
stream rapidly away leaving the chiral flelds in an unstable conflguration. During their 
subsequent evolution these flelds always "roll" towards the physical vacuum. This behaviour 
can be thought of as the inward collapse of the surface of the droplet as the quarks escape 
from it. Although such systems show no tendency to form DCC's, we found that coherent 
ampliflcation of any initial pion fluctuation could occur through a pion-laser-like mechanism. 
This is a consequence of the strong oscillations of the a fleld that can pump energy into 



oscillations of the pion fleld. Similar behaviour has also been seen in Refs. p3| , p^ . We have 
shown in I that this mechanism is a robust one: for example, it does not require a chiral 
phase transition in order to produce enhanced pion flelds. 

For cylindrical systems studied here, we flnd a competition between two mechanisms. 
The flrst consists of the transverse flow of the quarks and resulting collapse of the surface 
of the tube. This leads to behaviour that is similar to that seen in the spherical case. The 
second is the dilution of the quark density inside the tube as a result of its longitudinal 
expansion. The latter can lead to DCC formation, as one might expect from the results of 
Ref. IP]. However, the extent to which this mechanism operates is very sensitive to details 
of the initial conflguration. The formation time for any DCC relative to the time for the 
tube to collapse determines whether a DCC forms, and how long it can exist. The collapse 
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also changes the momentum distribution of the pions from that of the initial DCC, by 
coherently exciting pion modes with higher momenta. In addition we find that inclusion of 
explicit chiral symmetry breaking, with the strength needed to give the observed pion mass, 
significantly suppresses DCC formation. 

The paper is organized as follows. In Sec. II we present the linear sigma model that 
we use together with the classical transport and field equations. Since the basic approach 
has already been described in I, we concentrate on those features that are specific to the 
longitudinally expanding cylinder. Our results are described in Sec. Ill and we discuss their 
implications in Sec. IV. 



II. MODEL 



We work here with the linear sigma model ^7j, which provides a simple model for the 
physics associated with chiral symmetry. The model describes quarks interacting with a 
chiral four- vector of meson fields (a, tt). Since we study here configurations in which only 
one component of the pion field is nonvanishing, we keep only one pion field and, as in I, we 
further simplify the model by neglecting the isospin dependence of the quark-pion coupling. 
The Lagrangian we use is thus 



C = ^[i^-g{a + ^75)]^ + -{d^crd'^a + d^Tid'^n) - U{a, tt) 



(1) 



in which the meson fields of an 0(2) linear sigma model are coupled to two flavours and 
three colours of quark. 

The interactions among the meson fields are described by the potential U, which we take 
to be of the form 



4 



/ 2 I 2 2\2 

(a + TT — u ) 



(2) 



where /^r = 93 MeV is the pion decay constant. The "Mexican-hat" form of this potential 
leads to spontaneous breaking of the chiral symmetry through the nonzero vacuum expec- 
tation value of a, which corresponds to the quark condensate of the QCD vacuum. The 
parameter A is chosen to give the a meson a mass in the range 600-1000 MeV. In the phys- 
ical vacuum the quarks develop a mass Mg = gf^,. We consider values of g that correspond 
to quark masses in the range 300-500 MeV. 

For the situation considered here, it is convenient to work in terms of polar coordinates 
p and for the transverse position, proper time r = \Jt'^ — , and (space-time) rapidity 
7] = I ln[(t + z)/(t — 2;)]. For a cylindrically symmetric system undergoing a boost-invariant 
expansion, the fields depend only on p and r. In this case the equations for the meson fields 
take the form 





A2((T2(r,p) + 7r2(r,p) - z/^) + /5q(r, p)] a(r,p) + f^m\ 



(3) 




\ d ( d-K^ 
pdp V (9p , 



X^ia\T,p) + T:\T,p)-u')+g'S^{T,p) 7r(r,r). (4) 
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The source density in these equations, iSq(r, p), is proportional to the scalar density of quarks. 
Its detailed form is given below. 



The Vlasov equation is conveniently expressed in terms of momentum variables 
defined in the local rest frame of the matter. In this frame, which is specified by the unit 
four- vector u = {t,0,z)/T [^], the longitudinal momentum p\\ corresponding to a classical 
particle with three-momentum p moving in the presence of mean scalar and pseudoscalar 
fields is 

Pzt - Ez 

P\\ = ' (5) 

where 



E{x,p) = ^^p^ + M\x), (6) 

and 



M{x) = gsja^ix) + 'K^{x). (7) 

Similarly, the energy of such a particle in this frame is 

Et - p,z 

e = p-u = . (8j 

T 

The transverse components of the momentum, denoted p_L, are of course unchanged by 
the boost to the local rest frame. In terms of these variables, the Vlasov equation for the 
phase-space distribution of quarks /(r, 77, r_L,p||, p_|_) can be written 



d p\\ d I „ nil (9 , nN „ 



f{T,T],V±,Pl\,p^) = 0, 

(9) 



where we have assumed boost invariance, i.e. that M{x) is independent of 77. The antiquark 
distribution, denoted /(r, 77, r_L,p||, p_l), satisfies an equation of similar form. 

The Vlasov equation (^) describes freely streaming classical quarks and antiquarks. 
These particles have energies that are related to their three-momenta by Eq. and obey 
relativistic single-particle equations of motion, which, in terms of the coordinates and mo- 
menta defined above, take the form: 

rx(r) = HiW. (10) 
Pi{t) = -V,^(, (11) 

Mr) = ^. (12) 

re 

P||(r) = -^, (13) 
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where the overdots denote derivatives with respect to r and the particle's energy in the local 
rest frame is 



;(r, rx(r),p||(r), 



M2(r,r^(r))+p2(r) + p^(r). 



(14) 



The source term in the field equations, Sq{x), is given by the integral of [f{x,p) + 
f{x,p)]/ E{x,p) over all three-momenta [p^ , p!^] . Changing variables to p|| and p_L, we can 
write it in the form 



SJx) 



e(x,p||,px) 



(15) 



Rather than solving Eq. (^) directly as a partial differential equation in seven dimensions, 
we use the test-particle method [^,^]. In this approach, the smooth distributions / and / 
are approximated by a set of classical particles obeying the equations of motion (|10D-(|13[). 
For the numerical results presented here, we used 40 000 test quarks and antiquarks. 

The numerical techniques for solving the equations of motion (M), (H) and (Rl) are very 



similar to those applied to a soliton bag model in Ref. |2^. Details of the method for the 
case of the linear sigma model can be found in I. In the present case, the boost invariance of 
the system leads to a number of simplifications. The quark distributions are independent of 
7] and hence we do not need to consider the longitudinal motion of the particles, Eq. (P^. 
In addition, the longitudinal momenta satisfy Eq. (0) and so they simply scale like 1/r . 

Because of the longitudinal expansion, there is no finite, conserved energy for systems 
with this geometry. Nonetheless it is convenient to define an energy per unit rapidity in the 
local rest frame. This takes the form 



E 



a' + + \{y.,<yf + ^(Vr.vr)^ + f/(a, vr) 
y"rfp|| rf^pxe [/(r,r7,rx,P||,Px) + /(r, 77, rx,p||, p^) 



(16) 



where a constant has been added to t/(cr, vr) so that it vanishes in the vacuum. By making 
use of the equations of motion, one finds that the rate of change of E is given by 



dr J 12 



- i(Vr,a)2 - i(Vr,7r)2 - U{a, vr) 



+ Idpwd^p^— /(r,r/,rx,P||,p±) + /(r,r/,rx,P||,p±) 



(17) 



We have used this equation as a check on our numerical integration of the equations of 
motion. 



III. RESULTS 

In this section, we present the results of our simulations corresponding to various initial 
conditions of the system. As in I, the results shown are for our "standard" choice of pa- 
rameters, m^j = 1000 MeV, Mq = 300 MeV and, if chiral-symmetry breaking is included. 
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= 139 MeV. For these parameters, the temperature below which the phase with o" 7^ 
becomes the ground state is Tq ~ 235 MeV. 

The initial conditions are specified by three parameters: the radius tq of the tube, the 
temperature T of the quark plasma and the proper time tq at which the evolution according 
to Eqs. d), (^) and (|) starts. A chemical potential can also be introduced to allow for 
a nonzero initial baryon number, but this does not qualitatively change the behaviour of 
the system and so we present here only cases with /x = 0. We assume, as before, that the 
plasma inside the tube is initially uniform.]^ 

In order to study whether initial pionic fields can be amplified during the evolution, we 
add a small pionic perturbation to the initial configuration. As in our previous work, we 
have considered only uniform initial fiuctuations of the pion field since our aim is to study 
whether DCC formation is possible in these systems. A more realistic approach would be 
to take initial fiuctuations from a thermal distribution using a method similar to that in 
We would then expect to find the formation of regions of differently oriented DCC's, 



as seen in the work of Asakawa et al. |T^. Such studies will require integration of the full 
three-dimensional equations of motion and so will be much more computationally intensive. 

For many choices of initial conditions, and in particular for cases with relatively small 
initial tubes (rg < 5 fm), we find similar behaviour to that seen for the spherical droplets of 
quark plasma studied in I. In these cases, the rapid outward streaming of the quarks leaves 
the chiral fields in an unstable configuration. The surface of the tube collapses inwards, with 
the chiral fields "rolling" towards the physical vacuum. The a field then executes strong 
oscillations abouts its vacuum value. If a nonzero initial pion field is present, then this can 
be amplified by the laser-like mechanism also seen in I. 

In contrast, for large enough initial tubes, we find a rather different behaviour. This is 
particularly clear in the chiral limit. An example of this type is shown in Figs. |T]-^. These 
show the behaviour of a tube of initial radius tq = 6 fm and temperature T = 250 MeV 
at To = 1 fm/c. This initial proper time is the same as that used in Refs. [^,[TD|,|TT],|T^,|T3 



but is significantly smaller than the estimates of the freeze-out time in the work of Csernai, 
Mishustin and coworkers |^,18,20|. Note that the two examples for which we display results 



here are not necessarily the most realistic ones; rather they have been chosen since they 
most clearly illustrate the types of behaviour that are possible, and the conditions under 
which these can occur. 

Figs. |l| and § show the evolution of the chiral fields for this example of large tube with an 
initial temperature just above Tq. The behaviour of the central region of the tube is domi- 
nated, at least initially, by the longitudinal expansion. In particular the quark density there 
drops below its critical value before the quarks have a chance to start streaming outwards. 
At this point the chirally restored phase becomes unstable and any small fiuctuation from 
it can start to grow exponentially. In the present case, this happens at r ~ 2.5 fm/c, when 
the pion field at small radii can be seen to rise rapidly and then oscillate about a value close 
to 

This pion field is uniform across the central region of the tube. It is an example of a 
DCC: a region of misaligned vacuum. The collapse of the surface of the tube means that 



'For further details of the implementation of the initial conditions in this model, see Sec. IV of I. 
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the DCC does not persist in this form for proper times longer than tq/c. In Figs. |l] and 
one can see nonzero a fields appearing at successively smaller radii and the corresponding 
pion fields ceasing to oscillate in phase with that at the centre of the tube. From r ~ 7 fm/c 
onwards, the behaviour of the system resembles that of the ones studied in I, the sigma field 
oscillating violently until r ~ 12 fm/c and then settling down to its vacuum value. 

By Fourier analysing the pion field at successive times, we have also investigated the 
spectrum of the pion modes that are excited. In terms of the transverse Fourier transform 
of the field, 

7r(r, k_L) = j dh± exp(zk_L ■ r_L)7r(r, p), (18) 

and its proper-time derivative tt (r, k^) we define the corresponding intensity in momentum 
space (per unit rapidity): 

S^{t, k^) = I (|^(t, k^) p + oolmt, k^) |2) , (19) 

where Uk = \Jk']_ + m^. At large times, when the pion fields are sufficiently weak that they 
are well described by a linearised equation of motion, this is just the energy density of the 
pion field in momentum space. A convenient measure of the total strength of the pion field 
is provided by 

which, for weak fields, is equal to the total number of pions per unit rapidity. 

In Fig. ^we show the behaviour of Nt^{t). Initially this rises very rapidly and then shows 
strong oscillations. This is similar to what was seen in the studies of Rajagopal and Wilczek 
PJU. After r ~ 7 fm/c, these oscillations disappear and are replaced by a more gradual rise 
in Nt^{t). From the spectrum shown in Fig. |^, we see that the initial oscillations are due 
to the strong excitation of modes with the lowest frequencies and so are characteristic of a 
DCC. As the tube collapses, the nearly uniform DCC is destroyed and strength is removed 
from the lowest modes. Nonetheless significant coherent pion fields continue to be present 
although, unlike the original DCC, these oscillate in space and time. 

In Fig. ^ we plot the same spectrum as a function of pion momentum at various proper 
times. At r = 4 fm/c one sees the very large strength concentrated in modes with transverse 
momenta of less than 150 MeV. This pattern is characteristic of the initial, nearly uniform 
DCC. At later times, this strength falls off and modes with higher momenta become excited. 
In particular, significant strength builds up in modes with transverse momenta in the range 
300-500 MeV. There is also a smaller peak for momenta of 700-900 MeV Although the final 
pion fields are no longer dominated by the lowest momentum components, their amplitude 
still refiects the fact that a DCC was formed. For example, A^^ is enhanced by a factor of 
about 2000. This is far larger than the effects seen in the spherical systems studied in I, 
where the maximum enhancement factors were less than 100. 

When explicit chiral symmetry breaking is present, in the form of the final term in the 
potential (0), there is no phase transition but only a crossover between vacua with small 
and large values of the a field. For the physical value of m^r, the crossover is fairly rapid and 
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so some of the features of the phase transition do survive. Nonetheless exphcit symmetry 
breaking does lead to some qualitative differences in the evolution of the chiral fields after 
a quench. An example is shown in Figs. ^|-^ In this case the radius of the initial tube is 
again tq = 6 fm and tq = 1 fm. The initial temperature was taken to be T = 225 MeV, 
which lies within the crossover region for our standard parameter set. 

The explicit symmetry-breaking term in (|^) tilts the Mexican hat potential in the direc- 
tion of the true vacuum. Hence the system never evolves towards a maximally misaligned 
vacuum of the sort seen in the previous example. Despite this, significant amplification of 
any initial pion field can occur and so a DCC can still form. This is shown in Fig. |^ by 
the appearance of a more-or-less uniform pion field of ~ 20 MeV in the centre of the tube. 
This field starts to oscillate around the true vacuum, vr = 0. In an infinite system, like those 
studied in Refs. P,pO|, these oscillations would continue indefinitely. In the present case. 



they are terminated at r ~ 7 fm/c by the collapse of the tube. 

The fact that the fields inside the tube tend to roll towards the true vacuum, as a 
result of the tilted potential, also means that much less energy is released when the tube 
finally collapses. Indeed, in the example shown, the tube never completely collapses and a 
cylindrical region of cold quark matter is formed. Formation of such a region of quark matter 
was predicted by Csernai and Mishustin and was also seen in I for parameter sets with 
strong quark-meson couplings. This matter is created in a highly excited configuration, as 
can be seen in Fig. 6 from the strong oscillations of the a field from r ~ 10 fm/c onwards. As 
this matter settles down, it radiates a significant fraction of its energy in the form of pions, 
which is the reason for the continued growth of A^^r with proper time seen in Fig. |[ However, 
as discussed in I, we believe that the formation of such matter represents an artifact of the 
model and so this behaviour should not be taken too seriously. By r ~ 10 fm/c, before 
this behaviour sets in, A^^^ has been enhanced by a factor of about 200. Although this 
enhancement is significantly less than that in the previous example, it is still much larger 
than that found in the spherical systems studied in I. 



IV. DISCUSSION 

We have studied the evolution of systems of quarks coupled to chiral fields, starting from 
an initial tube of hot quark plasma undergoing a boost-invariant longitudinal expansion. 
Such a geometry is expected to be relevant to the central regions of ultra-relativistic heavy- 
ion collisions. We assume that a rapid quench occurs at some proper time tq, leaving the 
quarks out of thermal equilibrium with the chiral fields. The subsequent evolution of the 
system is treated in the classical approximation, the quarks being described by a relativistic 
Vlasov equation in the presence of a and pion fields which satisfy the mean-field equations 
of the linear sigma model. 

Unlike the spherical droplets studied previously in I, these cylindrical configurations can 
form DCC's of the sort seen in infinite systems following a quench Q. Similar behaviour is 



also seen in infinite systems undergoing longitudinal expansion [[4p|- [Tl|JT5|j20|] . However we 
find that the transverse flow of the quarks out of the tube tends to act against the formation 
of a DCC. This effect does not occur for infinite systems. It leads to an inward collapse of 
the surface of the tube, which is similar to what happens for spherical droplets. Although 
this collapse does not generate DCC's, it can enhance any initial pion fluctuation through 
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a pion-laser-like mechanism whereby strong oscillations of the a field pump energy into the 



pion fields |2|j2| . 

In the cylindrical case, the presence of two competing mechanisms leads to considerable 
sensitivity to the initial conditions. We find that significant DCC formation occurs only in 
large tubes (ro > 5 fm) where the longitudinal expansion rapidly reduces the quark density 
below the critical value needed to keep the system in the chirally symmetric phase. For a 
DCC to be produced, the initial density of quarks and antiquarks should be close to the 
critical value {i.e. the temperature should be close to Tq) and the proper time at freeze-out 
should be small compared to the radius of the tube. 

If the initial proper time Tq is taken to be ~ 1 fm/c, as in Refs. p,|T0|,|TT],|T^|T5[ and the 
illustrative examples shown here, then tubes with radii of ~ 5 fm or more have time to form 
significant DCC's. On the other hand, if tq is ~ 5 — 10 fm/c, as suggested by Csernai and 
Mishustin [^|18|,^, then it seems unlikely that DCC's can form in tubes of realistic radii. 



Clearly it is crucial to have better models for the early stages of ultra-relativistic heavy-ion 
collisions in order to determine the temperature and proper time at freeze-out. 

We also find that the subsequent evolution of such systems can significantly affect the 
transverse-momentum distribution of the emitted pions. Even if a system does form an 
initial DCC involving only very low-momentum modes of the pion field, the energy released 
by the oscillations of the a field as the tube collapses can lead to coherent excitation of 
pionic modes with higher transverse momenta. 

A further feature of our results is that explicit breaking of chiral symmetry has a sig- 
nificant effect on the evolution of these systems. In particular, by tilting the Mexican hat 
in the direction of the true vacuum, it reduces the degree of misalignment of the vacuum 
that can appear at the centre of the tube. A similar effect is also seen in infinite systems 
2^. As a result, the enhancement of the number of pions, as given by Eq. (^), is typically 



much smaller than that seen in the chiral limit. Nonetheless, even in the broken-symmetry 
example shown above, the pionic enhancement is much larger than that seen in the spherical 
systems studied in I, where no DCC was formed. 
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FIG. 1. Time evolution of the a field at various radii. The parameters of the model are 
m„ = 1000 MeV and Mq = 300 MeV, and the system is chirally symmetric {m.,^ = 0). The 
tube has initial radius ro = 6 fm and temperature T = 250 MeV at tq = 1 fm/c. 
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FIG. 2. Time evolution of the pion field at various radii. The model parameters and initial 
conditions are the same as in Fig. 1. 
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FIG. 3. Time evolution of Nt^. The model parameters and initial conditions are the same as in 
Fig. 1. 
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FIG. 4. Time evolution of the pion energy density in momentum space for the lowest ten modes 
of the field. The model parameters and initial conditions are the same as in Fig. 1. 
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FIG. 5. Momentum dependence of the pion energy density at various proper times. The model 
parameters and initial conditions are the same as in Fig. 1. 
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FIG. 6. Time evolution of the a field at various radii. The parameters of the model are 
= 1000 MeV, = 139 MeV and = 300 MeV. The tube has initial radius ro = 6 fm 
and temperature T = 225 MeV at tq = 1 fm/c. 
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FIG. 7. Time evolution of the pion field at various radii. The model parameters and initial 
conditions are the same as in Fig. 4. 
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FIG. 8. Time evolution of Nt^. The model parameters and initial conditions are the same as in 
Fig. 4. 
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